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ABSTRACT 

We present high-resolution radio images at 1.4 GHz of two Galactic 
supernova remnants (SNRs), G003.8-00.3 (formerly G003. 7-00.2) and G350.0- 
02.0 (formerly G350.0-01.8). Although the two objects are very different in 
appearance, in both cases the radio emission shows a clear bilateral (or "barrel" ) 
morphology for which the axis is parallel to the Galactic Plane. 

The majority of Galactic SNRs have now been observed at high resolution, 
and one can define a clear bilateral subset of the population. We consider a 
sample of 17 such SNRs, and find a highly significant tendency for the bilateral 
axes of these SNRs to be aligned with the Galactic Plane. We interpret this as 
indicating that "extrinsic" effects dominate the morphology of such remnants. 
Specifically, we argue that the Galactic magnetic field causes these SNRs to 
appear bilateral, either directly, in the form of magnetic field compression 
and/or quasi-perpendicular acceleration of electrons in the supernova shock, 
or more likely indirectly, by pre-processing the interstellar medium to produce 
density stratifications extended along the Plane. 



Subject headings: Galaxy: structure — ISM: magnetic fields, structure — shock 
waves — supernova remnants: individual (G003.8-00.3, G350. 0-02.0) 
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1. Introduction 

High resolution observations of radio supernova remnants (SNRs) have revealed 
exceedingly complex structures. The appearance of a SNR can tell us much about the 
supernova (SN) explosion itself, the nature of any central compact object, the circumstellar 
medium (CSM) of the progenitor, the surrounding interstellar medium (ISM), and the 
properties of the ambient magnetic field. Indeed the difficulty in interpreting SNR 
morphologies has always been in disentangling the competing effects and in determining 
which ones dominate during the different evolutionary stages. 

Perhaps the longest running debate concerning SNR morphologies has been over the 
classification referred to as bilateral, bipolar, axially symmetric, or "barrel" SNRs. These 
remnants are characterised by a clear axis of symmetry, low levels of emission along this 
axis, and two bright limbs on either side (e.g. pardner fc Milne 1965| ; [Whiteoak fc Gardner 



19681 [Roger fc Costain 197q ; |Kesteven fc Caswell 1987| , hereafter KC87; [Roger et al. 198^ 



[Caswell et al. 1992|) . It is generally accepted that bilateral SNRs represent an underlying 
cylindrical symmetry, the emitting regions corresponding to the curved walls of the cylinder 
(KC87). A bilateral appearance is then produced when the "barrel" axis is approximately 
in the plane of the sky, the two bright fianks running parallel to this axis. 

Many mechanisms for generating a bilateral appearance have been suggested (see 
Bisnovatyi-Kogan fc Silich 1995[ for a review). However it is not clear which of these models. 



if any, actually apply. Indeed, given the diversity in the environments and ages of Galactic 
SNRs, several of the mechanisms proposed may be valid. We emphasise that observations 
demonstrate that many SNRs are not bilateral. To some extent, this can be accounted for 
by viewing angle (see KC87 and Section p.4[ below), but it is important to keep in mind 
that whatever mechanism is invoked, the morphology may be dominated by small scale 
effects such as the inhomogeneity of the local environment. 

We now briefiy outline previously proposed explanations for bilateral SNRs. 
Explanations can be divided up into two categories: "extrinsic", relating to the ambient 
ISM and magnetic field, and "intrinsic", relating to the progenitor, its CSM, the SN 
explosion itself, and any compact stellar remnant. 



1.1. Extrinsic explanations 

The simplest model to explain a bilateral appearance is that it is the result of the 
density structure of the surrounding ISM. Simulations show that when a remnant expands 
inside an elongated, or tube-like, cavity, it quickly takes the shape of that cavity, and a 
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bilateral appearance can be produced ( [Bisnovatyi-Kogan, Lozinskaya, fc Silich 1990 ). 

It has long been realised that the geometry of the ambient magnetic field can also 
affect the appearance of a SNR. We collectively refer to models relying on the ambient field 
as "magnetic" . In a uniform magnetic field, the shock front produced by an expanding SNR 
will preferentially compresses the field where the shock normal is perpendicular to the field 
lines. Thus, in a well-ordered ambient field with a significant component in the plane of 
the sky, enhanced emission along two limbs is generated, producing a bilateral appearance 
( [van der Laan 19"6^ ; [Whiteoak fc Gardner 19681) . 



Argument is divided as to whether this mechanism can produce the observed ratio 
of maximum to minimum intensity around the circumference of SNRs which are in the 
adiabatic phase of evolution ([Fulbright fc Reynolds 1990 , hereafter FR90; P^atkiewicz 



Axford, fc McKenzie 199^ . Thus, an extension of this model has been proposed, in which 
the efficiency of particle acceleration by the SN shock depends on the angle between the 
shock normal and the orientation of the ambient magnetic field ( Jokipii 1987] ; [Roger et al. 



1988[ ; FR90). FR90 argue that the bilateral appearance in the adiabatic phase can best 



be explained by "quasi-perpendicular" acceleration, in which acceleration is most efficient 
(and synchrotron emission strongest) where the shock normal and the field are at right 
angles. "Quasi-parallel" models (where the shock normal and the field are aligned) can 
also produce a bilateral appearance at some orientations, but at other orientations produce 
morphologies which are never observed (FR90, but see [Zhang, Wang, fc Chen 1996[) . 



1.2. Intrinsic explanations 

There are a number of symmetry axes associated with the progenitor system and the 
explosion. Several mechanisms have been suggested by which these symmetries can produce 
a bilateral appearance. Models include a toroidal distribution of ejecta (Podenheimer[ 
Woosley 1983[ ; KC87), the effect of a high velocity progenitor ( [Rozyczka et al. 1993[ ; 



Prighenti fc D'Ercole 199^ ), the distribution of mass loss and magnetic field in the CSM 
produced by the progenitor ( [Manchester 1987| ; [Igumenshchev, Tutokov, fc Shustov 1992[ ; 



[Storey et al. 1992[ ; [Zhang, Wang, fc Chen 1996| ), and the infiuence of outfiows from a 
central compact object ( [Manchester fc Durdin 1983[ ; [Willingale et al. 1996[ ). 
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1.3. A statistical approach 

Understanding the morphology of an individual SNR can be difficult, often requiring a 
concerted multi-wavelength campaign. Only a few bilateral SNRs can be considered to be 
well-studied (most notably G296. 5+10.0 and G327. 6-1-14. 6) and even these objects provoke 
controversy (e.g. compare Roger et al. 1988 with Storey et al. 1992). An alternative 
approach is to consider the common properties of bilateral SNRs on a statistical basis. 

There is reasonable evidence that at least for low latitudes, the Galactic magnetic field 
is well-ordered and generally runs parallel to the Galactic Plane when projected onto the 
sky ( llVlathewson fc l-brd 1970| ; |lVlanchesterT97^ ; [Ellis fc Axon 1978| ; [Sofue fc l-ujimoto] 



[Rand fc Kulkarni 1989| ; [Reid fc Silverstein 1990| ; see also reviews by [Sofue, I'ujimoto 



Wielebinski 19"86| and [Beck et al. 1996| ) . If a bilateral appearance is caused by a magnetic 



model, one might expect the structure of the Galactic magnetic field to manifest itself as 
a statistical alignment of the symmetry axes of bilateral SNRs with the Galactic Plane 
( [Shaver 1969|) . 



Given that such a result would be reasonably conclusive, several researchers have 
attempted to determine whether such an alignment exists. By creating an average 
circumferential profile from SNRs for which appropriate data was available at the time. 
Shaver (1969) found that SNRs appeared to have minima in their emission where the shell 
is intersected by a line of constant Galactic latitude, and argued that this was evidence for 
magnetic models. However, the statistical significance of the result was low, and the limited 
sample available was from observations of poor resolution and sensitivity. 

In more recent years, the statistics of SNR samples have steadily improved, and other 
researchers have carried out various analyses investigating alignment between the axes of 
SNRs (defined in various ways) and the Galactic Plane ( Shaver 1982| ; [IVlanchester 1987 ; 



Leckband, Spangler, fc Cairns 1989| ; [Whiteoak fc Green 1996| , hereafter WG96). In all cases. 



no alignment was found, leaving the authors to conclude that either the magnetic field is 
not well-aligned with the Plane at such scales, or that the magnetic model is not valid. 

However, some of these studies implicitly assume that all SNRs have a bilateral 
symmetry, and then assign an axis to every remnant in their sample. Others restrict their 
sample to SNRs with a symmetry axis, but include remnants where this symmetry is not 
necessarily bilateral. If only a small subset of remnants are truly bilateral, then these 
analyses will never show any alignment because any real effect will be swamped by the 
large number of randomly oriented axes corresponding to non-bilateral SNRs. In addition, 
the low resolution of many of the images available at the time caused a poor determination 
of the axis of symmetry. 
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There are now over 200 SNRs known in the Galaxy ( Green 1996 ), and an increasing 
number of objects have been observed at high resolution. For almost the entire Galactic 
population, observations are now of sufficient quality that one can determine which SNRs 
have a clear bilateral appearance (defined in some manner), and restrict analysis to this 
sample. 

In Section ^, we describe high-resolution 1.4 GHz observations and analysis procedures 
for two Southern SNRs, G003.8-00.3 (formerly G003. 7-00.2) and G350.0-02.0 (formerly 
G350.0-01.8). We present images of these remnants in Section ^, and argue that these 
SNRs are members of the bilateral class with their axes aligned along the Galactic Plane. 
The properties of these two SNRs are then briefiy discussed in Section ^. In Section ^ 
we consider the Galactic population of bilateral SNRs, and demonstrate a clear tendency 
for the bilateral axis to align with the Plane. It is argued in Section ^ that the magnetic 
model can explain this alignment, but that there are difficulties accounting for the detailed 
morphologies of bilateral SNRs. An alternative model is proposed, in which the ISM is 
pre-processed by the ambient magnetic field to form structures running parallel to the 
Galactic Plane. SNRs then interact with this material, producing a bilateral appearance 
with the observed alignment. 



2. Observations and data reduction 

Interferometric observations were made with the Very Large Array (VLA) of the 
National Radio Astronomy Observatory^, an aperture synthesis telescope consisting of 27 
antennae in a "Y" -shaped array, located near Socorro, New Mexico ([Napier, 'I'hompson, 



fc Ekers 1983| ). The dates and lengths of the observations are given in Table |1]. Two 



bands were observed simultaneously, with centre frequencies of 1.385 GHz and 1.465 GHz, 
bandwidths of 50 MHz, in all four Stokes parameters. Because of its large angular extent, 
SNR G350. 0-02.0 was observed in a mosaic of four pointings (see Table |^), with a cycle 
time of 24 minutes. Amplitudes were calibrated by assuming fiux densities at 1.385 and 
1.465 GHz respectively of 14.94 and 14.55 Jy for J1331+305 (3C286) and 16.32 and 15.62 
Jy for J0137-I-331 (3C48), where 1 Jy = 10~^^ W m~^ Hz~^. Phases were calibrated using 
the source J1751-253. 

The VLA data were edited and then calibrated in AIPS using standard techniques 
( IGreisen 1996|) . The visibility data were then transferred to the Miriad package (|Sault,| 



Jrhe National Radio Astronomy Observatory is a facility of the National Science Foundation operated 



under cooperative agreement by Associated Universities, Inc. 
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Teuben, & Wright 1995), and total intensity images formed using uniform weighting and 
a cell size of 4". For the mosaiced observations of G350. 0-02.0, a joint image was formed 
from all four pointings ( |(Jornwell, Holdaway, fc Uson 1993| ; |Sault, Staveley-Smith, fc Brouw 



1996 ), which was then corrected for the mean primary beam response of the VLA antennae. 

The smallest projected spacing between antennae for the observations in Table |l| is 
0.12 kA , which corresponds to an angular scale of ~ 30'. Thus the VLA observations are 
not sensitive to emission on scales larger than this. These short spacings were obtained by 
making single-dish observations (in total intensity only) of both SNRs using the Parkes 64m 
radio telescope[[ Amplitudes were calibrated using the source PKS B1934-638, assuming a 
flux density of 14.90 Jy. Except where noted in Table 0, observational techniques were the 
same as used in the 2.4 GHz survey of Duncan et al. (1995). Scans were made in both / and 
6, of regions surrounding both SNRs, at a scan rate of 6° min^^. The Parkes observations 
were reduced in the manner described by Duncan et al. (1995), and the resultant images 
regridded and reprojected to match the higher resolution data. 

Data from the two telescopes were then combined by adding together the dirty images 
from each instrument. The flux density scale of the single dish data was preserved by 
weighting them by the inverse ratio of the beam areas ( Ye, I'urtle, fc Kennicutt Jr. 199T| ; 



Stewart et al. 19"93|) . Point spread functions were produced similarly, and the images then 



deconvolved using the maximum entropy algorithm ( |Gull fc Daniell 1978|) . For G350. 0-02.0, 



all pointings were deconvolved simultaneously using the Miriad task MOSMEM (jSault 



Staveley-Smith, fc Brouw 1996|) . The resulting models were then smoothed using a Gaussian 



restoring beam to produce the final images. 

Because of the proximity of G003.8-00.3 to both the Galactic Plane and the Galactic 
Centre, single-dish observations of this source are dominated by diffuse emission from 
confusing sources and the Galactic background. The composite image has flux filling the 
entire field of view, and it is thus not possible to successfully deconvolve it. The image of 
G003.7-00.2 as presented in Section || therefore does not include the Parkes measurements. 
Since the angular extent of this SNR is significantly smaller than the largest spatial scale 
sampled by the VLA, we still expect to recover almost all of the flux density from this 
object. 

Images of linearly polarized emission were made from the interferometric data in Stokes 
Q, U and V. (Although no real emission is expected in circular polarization, the V image 
provides a useful check of quality.) The Q and U images were then deconvolved using the 



*The Parkes radio telescope forms part of the Australia Telescope, which is funded by the Commonwealth 
of Australia for operation as a National Facility managed by CSIRO. 



- 7- 



CLEAN algorithm ( |Clark 1980 ), and a linear polarization image, L, produced which was 
then corrected for non-Gaussian noise statistics ( [Killeen, Bicknell, fc Ekers 1986| ) using the 
Miriad task IMPOL. 



3. Results 

Total intensity images of G003.8-00.3 and G350. 0-02.0 are shown in Figures |1] and ^ 
respectively. These images are of higher resolution, and in the latter case, of greater extent 
than in previous images. We are thus able to derive improved positions for the centres of 
these two SNRs, resulting in a change of name for both remnants. 

Because the VLA is non-coplanar, the approximation of a flat sky used in standard 
synthesis imaging is invalid at large distances from the phase centre ( Perley 1989| ). As a 
result, some point sources in Figure ^have a triangular appearance. 

Information about each image and its derived parameters are provided in Table |^. Flux 
densities and their uncertainties were determined by integrating within multiple polygons 
enclosing the emission from each SNR. The RMS noise in each image was measured 
by similarly integrating over a source-free region. In the case of G003.8-00.3, the given 
centre and diameter are those of a circle best fitting the outer perimeter of the emission. 
G350. 0-02.0 is clearly non-circular, and so the centre is defined as that point equidistant 
from the boundaries of the SNR both along and at right angles to the bilateral axis. The 
diameter given is the maximum extent in these two directions. 



3.1. G003.8-00.3 

SNR G003.8-00.3 (then G003. 7-00.2) was discovered in the Galactic Centre Survey 
carried out with the Molonglo Observatory Synthesis Telescope (MOST) at 0.843 GHz 
( |Gray 19"^ ) and was classified as a SNR on the basis of its morphology. Using flux 
densities from Gray (1994) and from Table ^, we obtain a spectral index for G003.8-00.3 of 
a = —0.65 ± 0.05 {S^, oc z/"), confirming the non-thermal nature of this source. 

The image of G003.8-00.3 (Figure P shows it to be a classic bilateral remnant, with a 
clear axis of symmetry defined by the two bright limbs, and a low level of emission along 
this axis (cf. KC87). The emitting regions appear to be composed of multiple filaments 
extended along the symmetry axis: on the north-west rim of the SNR, these filaments take 
the form of roughly parallel arcs, while in the south-east they overlap each other. The 
original 0.843 GHz image of this remnant seemed to contain multiple overlapping rings 
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( CJray 1991 ). However this appearance is most probably due to the ahgnment of various 



filaments when seen with the lower resolution of the MOST — no such features are evident 
in our image. 

ifj is defined to be the acute angle between the symmetry axis of the SNR and the 
Galactic Plane. We fit the symmetry axis by eye, and estimate the uncertainty by fitting a 
set of possible axes to the remnant. This gives a value ip = 6° ±2°, indicating close alignment 
between the symmetry axis and the Plane. More involved methods for determining ip (e.g. 
^torey et al. 199"^ ) were also implemented, giving similar results. For the accuracy required 
for ip in this work, the more rigorous determinations are not warranted. Furthermore, for 



several sources considered in Section p.3|, images were not available in electronic form. 



Weak linearly polarized emission fills the entire image, and it is not possible to 
determine how much emission is specifically associated with the remnant. We put an upper 
limit of 0.05 Jy on polarized emission from G003.8-00.3, corresponding to a fractional 
polarization of <3%. 



3.2. G350.0-02.0 

SNR G350.0-02.0 (then G350.0-01.8) was first identified from a 408 MHz survey of the 
Galactic Plane carried out with the Molonglo Cross ( [Green 1972| ). At the comparatively 



low resolution and sensitivity of the time, only the bright north-western arc was clearly 
delineated. For the next 20 years, this remnant received little attention: Clark et al. (1975) 
suggested it was "part of an old large SNR", Caswell (1977) noted it as an example of a 
SNR which was brighter towards the Galactic Plane, while KC87 described it as "difficult 
to classify". More recent observations at 0.843 GHz (MOST; [Burn fc Bush 1994| ) and at 2.4 
GHz (Parkes 64m; Puncan et al. 1995|) have demonstrated that this SNR is of far greater 



extent than previously realised. Higher resolution VLA observations are shown in Figure ^ 

Using a 2.4 GHz flux density of 18 ± 1 Jy (obtained from the single-dish measurements 
of [Duncan et al. 1995| ), and the 1.4 GHz value given in Table ^, we derive a spectral index 



for G350. 0-02.0 of a = -0.4 ± 0.1. Note that the flux density values quoted for this SNR 
by Clark et al. (1975) correspond only to the bright north-west component. 

The image of this SNR in Figure |^ shows a far more complicated morphology than 
previously thought. It may even correspond to multiple remnants, a possibility discussed 
below in Section [4.2| . The brightest region of G350. 0-02.0 (corresponding to the previously 



identifled source G350.0-01.8, and labelled as A in Figure 0) consists of several overlapping 
fllaments, with the suggestion of complex, unresolved, structure. This region, which extends 
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over ~ 40', forms a roughly circular arc with radius of curvature ~ 55'. In the remnant's 
interior is region B, separated from region A by a zone of diffuse emission of width ~ 8'. 
Region B consists of a series of wisps, considerably fainter than the emission in region A, but 
with a similar orientation and curvature. A third region, C, is separated from region B by 
another zone of diffuse emission of width ~ 6'. Region C consists of still more overlapping 
filaments, of opposite curvature to the rest of the SNR, which complete a distorted shell 
of emission. The outer rim of region C forms a circular arc of radius ~ 14', considerably 
different from the radius of curvature for region A. 

Although G350. 0-02.0 is clearly very different to G003.8-00.3, it also has a bilateral 
appearance. The circumference of the remnant is bounded on two sides by regions of 
significant emission (namely A and C), and at perpendicular position angles the emission is 
considerably weaker. As for G003.8-00.3, a bilateral axis can be defined. This axis is shown 
in Figure drawn perpendicular to a line joining the centres of curvature of regions A and 
C. For this axis, i/j = 5° ±2°, again a high degree of alignment with the Plane. 

Clumps of linear polarization are detected throughout the interior of G350. 0-02.0. The 
total polarized intensity is 0.8 Jy, implying a fractional intensity of 3.6%. 



Discussion 



Because of the lack of information at other wavelengths and the proximity of both 
remnants to / = 0°, the methods typically used to estimate distances to SNRs cannot 
be applied (see |Green 19B^ ). The S — D relation (e.g. [Gaswell fc Lerche 1979| ; pduang 



fe 'I'haddeus 1985| ) has in the past been used to estimate distances to SNRs, but the 
uncertainties are large ( preen 1984| ; [Berkhuiisen 1986| ). Thus we conclude that the age of 
and distance to these SNRs remains unknown at the present time. 



4.1. G003.8-00.3 

FR90 have quantified the contrast between the brightest and faintest regions of a 
bilateral SNR by using an "azimuthal intensity ratio", A. This is defined as the ratio 
between the peak of emission around the SNR shell, and the minimum value around the 
shell at the radius corresponding to this peak. The filamentary nature of this SNR makes 
it difficult to compute A in this manner, but we have made a commensurate calculation as 
follows. The shell of emission from the SNR lies within an annulus of inner radius 4'.2 and 
outer radius 6'.7, centred on the coordinates for G003.8-00.3 as listed in Table |[ We have 
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divided this annulus into 360 azimuthal bins, and computed the average flux density in each 
bin. A peak of approximately 1.1 mJy beam~^ occurs at position angles (measured north 
through east) of 100° and 260°. The minimum level of emission occurs at PA ^ 40°, where 
no emission can be detected above the 3cr noise level. This corresponds to an intensity ratio 
of ^ > 6, at a resolution of ~70 beams per diameter. For a remnant in the adiabatic phase, 
this is consistent with magnetic models for a bilateral appearance, but does not distinguish 
between field compression alone and quasi-perpendicular acceleration (FR90). 



4.2. G350.0-02.0 

One interpretation for the three spatially and morphologically distinct regions of 
emission in G350. 0-02.0 is that they are two or three separate SNRs, as has been postulated 
in objects such as G053.6-02.2 (3C400.2; Pubner et ai 1994| ) and 0547-697 (DEM L 316; 



Williams et al. 1997| ). Certainly regions B and C are of similar surface brightness, which 



could be interpreted as a complete shell, separate from region A. 

If we consider a region on the sky defined by 330° = —30° < / < 30°, 1.5° < |6| < 3° 
(i.e. longitudes and latitudes representative of G350. 0-02.0), we find nine SNRs in the 
catalogue of Green (1996). Given this sky density, the probability of finding two randomly 
distributed remnants less than 30' apart (as required here) is 14%. However, this is 
probably a lower limit, since the small number of SNRs detected in this region may be 
due to the incomplete sky coverage of existing surveys. Also, SNRs are not necessarily 
randomly distributed, but at least for those formed from massive stars probably represent 
the clustering of their progenitors. Thus there is no convincing statistical argument that 
G350.0-02.0 is a single SNR. 

Nevertheless, the centres of curvature for all three regions lie in a straight line, which 
bisects the three arcs of emission. The fact that regions A and C are of opposite concavity 
also suggests that all the emission represents common expansion from a single event. Diffuse 
emission from the SNR (contributed primarily by the single dish observations) delineates 
a single plateau of emission bounded by regions A and C. We therefore consider the most 
likely interpretation to be that the emission in this region is from a single object. 

It is possible that the non-circular appearance could be due to the influence of a pulsar 
or other compact object. An analogy could be drawn with SNR G005.4-01.2, which has 
a bright, flattened ridge of emission along one side ( [Caswell et al. 1987|; [Frail, Kassim, fc 



[Weiler 1994[) , resembling region A of G350. 0-02.0. It has been proposed that the distorted 
morphology of this SNR is a result of the associated pulsar PSR B1757-24 overtaking the 
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expanding parent shell ( It* rail fc Kulkarni 1991 ). While a fiat-spectrum region extending 



from the pulsar back towards the radio shell supports such a model for SNR G005.4-1.2 
( [brail fc Kulkarni 1991| ; [brail, Kassim, fc Weiler 1994| ), no similar component is observed 



in G350. 0-02.0, nor has any pulsar been detected ( Manchester, D'Amico, fc 'I'uohy 1985[ 



[Kaspi et al. 19961) . Thus we consider an interaction with a compact object an unlikely 



explanation for the observed morphology. Given that any associated pulsar should be 
located near the SNR's presumed centre of expansion ( Paensler fc Johnston 1995| ), the area 
bounded by regions B and C might be a more fruitful target for future pulsar searches. 



A 1.4 GHz image of SNR G166.0+04.3 (VRO 42.05.01; [Landecker et at. 1982[ ; [Fineault 



Landecker, fc Routledge 1987[ ) is shown in Figure ^, revealing a striking resemblance 



between this SNR and G350. 0-02.0. G166.0+04.3 is also composed of three separate 
regions, whose relative shape, size and orientation are similar to those in G350. 0-02.0. Both 
SNRs are bounded on one side by a bright arc with a large radius of curvature which trails 
into the noise (region A), contain a fainter strip of emission parallel to A and separated 
from it (region B), and a series of filaments with a distinctly smaller radius of curvature and 
of opposite concavity to the rest of the emission, completing the outline of a shell (region 
C). Both SNRs are similarly aligned to the Galactic Plane: when an axis of symmetry is 
defined for G166.0+4.3, one finds that ip = 2° ± 2°. 

There are notable differences between the SNRs, however. Firstly, in G350. 0-02.0, 
regions B and C are of similar brightness but are much fainter than A, while in G166.0+04.3, 
the region corresponding to C (the "shell" of [Landecker et al. 1982| ) is of comparable 
brightness to that corresponding to region A (the "wing"). Secondly, the central region of 
G166.0+04.3 consists of linear features, while region B appears to be curved. Finally, while 
in G350. 0-02.0 regions B and C are separated by a region of faint emission, the analogues 
of these two regions in G166.0+04.3 actually overlap. 

Pineault et al. (1987) explain the morphology of G166.0+04.3 in terms of the 
interaction of an expanding remnant with a slab of pre-existing low density material. The 
SN explosion is believed to have occurred at or near the centre of curvature of the "shell" 
component, which represents the expansion of that part of the SN shock into a uniform 
medium. The faint linear features in the remnant's centre mark the re-energising of the 
slab's surface as a result of the shock breaking out into it, and the "wing" component 
represents the shock once again encountering dense material on the other side. 

The similarity between the two remnants suggests a similar explanation for them as 
well. In this model, region C of G350. 0-02.0 is the unperturbed component of the original 
SNR, which evolves independently of those parts of the remnant which encounter a more 
complicated structure ([Falle &: Garlick 1982[ ; [Tenorio-Tagle, Bodenheimer, fc Yorke 1985[) . 
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Indeed, region C resembles SNR G119.5+10.2 (CTA 1; [Pineault et al. 1993| ), a partial shell 
for which there is evidence that a breakout has occurred but for which the remainder of the 
SN shock is not visible. Region B of G350. 0-02.0 represents the boundary of a low density 
slab. That it is concave and not linear and does not overlap with region C, suggests that 
the details of the 3D geometry differ from that in the case of G166.0+04.3. Region A is that 
component of the shock which has traversed the slab and continues to propagate beyond it. 

5. The Galactic population of bilateral SNRs 

5.1. Selection criteria 

Although previous studies have argued that no alignment exists between the axes of 
bilateral SNRs and the Galactic Plane, it is interesting that G003.8-00.3, G350. 0-02.0 and 
G166.0+04.3 should all have such low values of iIj. This similarity prompted an examination 
of other bilateral SNRs. 

We have therefore reviewed the Galactic SNR catalogue of Green (1996), which 
currently contains 215 remnants. The criteria for a SNR to be classed as bilateral are as 
follows: 

• It must be of the shell or composite class. 

• The highest resolution image available must have a minimum of 10 beams across its 
diameter. 

• The SNR must have clear minima in emission separated by position angles of 
180° ±30° relative to the assumed centre of the SNR. The presence of opposed minima 
is the main criterion for classification, as it is these minima which define the bilateral 
axis. 

• The SNR must have well-defined maxima, similarly opposed, and at approximately 
perpendicular position angles to the minima. 

• A clear bilateral axis should be identifiable, passing through the two minima and 
through the centre of the SNR. Objects with a large uncertainty 20°) in the 
position angle of the axis are excluded — this requirement is critical in defining the 
objects in the sample. 

We do not quantify the azimuthal intensity ratio. A, of FR90 as a threshold. Such 
a measurement is both difficult and misleading, because of the various presentations of 
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images in the literature, and of the differing nature of the observations involved. As FR90 
have pointed out, the value of A depends strongly on the resolution of an image, so that 
the appearance of bilaterality will increase with resolution. Also, the lack of large-scale 
emission present in some interferometric images can cause a further over-estimation of this 
parameter. Thus using this ratio as a threshold will not result in a consistent criterion. 



5.2. Non-bilateral SNRs 

There are 198 SNRs in Green's (1996) SNR catalogue which do not meet our criteria. 
This is in strong contrast to KC87, who found that 63% of all SNRs had some element 
of "barrel" structure. This primarily reflects the strictness of our criteria. For example, 
consider G315.4-02.3 (ROW 86; Figure D and G332.4-00.4 (RCW 103; Figure |), both 
of which are described by KC87 as "well developed barrels" . The former has no opposed 
maxima or minima. The latter is a marginal case (Dickel et al. 1996 describe it is an 
"incipient barrel"), but lacks the distinctive minima and symmetry properties of the best 
examples of the class. 

Improved observations over the last ten years have allowed a more precise classification 
of many remnants. For the MOST, substantial improvements to both hardware (Amy & 
Large 1990, 1992) and data reduction techniques ( Pram fc Ye 1995| ) since the observations of 



KC87 were made have significantly reduced the effect of artifacts, increased the sensitivity, 
and generally improved image quality. 



5.3. Bilateral SNRs 

There remain 17 remnants in the catalogue which can be classed as clearly bilateral 
from our criteria. A list of these bilateral SNRs is given in Table ^ For each remnant a 
bilateral axis is defined by eye, as for the SNRs in Section ^. This axis passes through the 
two minima, corresponding closely to an axis of mirror symmetry. Some examples of SNRs 
which we have classed as bilateral are shown in Figures |^ through |TT|, with the bilateral 
axis of each remnant shown. In each case, we have measured the acute angle, ip, between 
this axis and the Galactic Plane. Values of i/j for each SNR are shown in Table 

The reader is encouraged to examine these (and other) SNRs in order to verify that 
the classifications made are reasonable. We should emphasise that if there is any doubt, a 
SNR is excluded from our class. If a preferential alignment exists, it will remain apparent 
if some bilateral SNRs are excluded from the sample (provided that the statistics are still 
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meaningful), but will be masked if too many remnants which do not satisfy the criteria are 
included. 



5.4. Orientation with respect to the Galactic Plane 

A histogram of the distribution in ifj from Table ^ is shown in Figure |12|. This 
distribution has a median of ipm = 12°l2°, where the uncertainty corresponds to the 
maximum shift in median should two objects be discarded from the sample. This result 
suggests an alignment of bilateral SNRs with the Galactic Plane: the probability of this 
distribution occurring by chance is now discussed. 

Assume cylindrical symmetry for bilateral SNRs, and represent the three-dimensional 
axis of such a remnant by a vector in Cartesian space, with spherical coordinates 9, (p as 



shown in Figure 0. The observer views the system along the r axis, and the / and b axes 
represent increasing Galactic longitude and latitude respectively in the plane of the sky. 
Note that the observed axis is a projection onto the / — b plane of the true axis, and that 
a remnant will only have a bilateral appearance provided that 6 > 6c, where 9c depends 
on the details of the geometry (KC87). We only consider values in the range < 90°, 
since the other hemisphere is equivalent. Randomly oriented vectors will have a uniform 
distribution in (p and hence in ip. Thus all values of tp are equally likely, independent of Re- 
using the binomial theorem, we can calculate that the probability of a randomly 
oriented sample having the observed median is CTlo^i x 10~^. We argue in Section |6.^ 
that G296. 5+10.0 and G327.6+14.6 may be distinct from the other objects in Table ^ as a 
consequence of their large distance from the Plane. These two remnants have the highest 
values of ip, and if they are excluded from the sample, the probability is even more tightly 
constrained. 

This calculation shows that the probability of the distribution being randomly oriented 
is low. We therefore consider our data to demonstrate strong evidence that bilateral SNRs 
align with the Galactic Plane. 

Given that no correlation between the Galactic Plane and the rotational axes of SNR 
progenitors has been observed ( ^lettebak 1949| ; [Huang fc Struve 1954|) , intrinsic models 



for the bilateral appearance (as listed in Section |1.2|) are not consistent with the observed 
distribution of ip- These models are probably better suited to explain the appearance of 
other classes of remnant morphology. For example, a toroidal explosion may be appropriate 
for G292.0+01.8 ( 'I'uohy, Clark, fc Burton 1982 ), a high velocity progenitor could describe 



G357. 7-00.1 ( [Shaver et al. 1985| ; [Becker fc Helfand 1985| ), outflows from a central source are 
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appropriate for G039. 7-02.0 (W50 and SS433; [Elston fc Baum 1987|) , while the properties 
of the GSM perhaps only cause a bilateral appearance at an early age, such as is observed 
in the radio remnant of SN 1987A (|Gaensler et al. 19971). 



6. Models for bilateral SNRs 

6.1. A magnetic model 

As discussed in Section |T]^, an alignment of bilateral SNRs can be interpreted as 
evidence for magnetic models ( ^haver 1969|; [Manchester 1987] ). In order for such mechanisms 



to be effective, the surrounding magnetic field must be well-ordered and the ISM into which 
the SNR is expanding must be reasonably homogeneous. The small number of bilateral 
remnants may indicate that this condition rarely occurs, as will be discussed further below. 

There is some support for the magnetic model from polarization measurements. 
Specifically, in the field compression scenario ( [Whiteoak fc Gardner 19BE| ), one expects 
the magnetic field direction in the emitting regions to be oriented tangentially. Such a 
phenomenon is observed in the bilateral SNRs G296. 5+10.0 ([Milne fc Haynes 1994|) and 



G093.3+06.9 ( [Lalitha et al. 1984[ ). A tangential magnetic field need not always be observed. 



however, since in younger SNRs other factors may cause the magnetic field to be oriented 
radially ( [Reynolds fc Gilmore 199"3[ ; [Jun fc Norman 1996[) or be generally disordered (FR90; 



[Mottett fc Reynolds 1991 



One argument against magnetic models is that the typical energy density in the 
ambient magnetic field is many orders of magnitude less than the kinetic energy of a typical 
SNR shock ([Manchester 1987t Pisnovatyi-Kogan, Lozinskaya, fc Silich 1990| ). However 



these arguments only pertain to whether the shape of a remnant can be influenced by the 
surrounding field. The brightness distribution can be significantly affected by the ambient 
field, even if the field is dynamically unimportant (FR90). 



6.2. An alternative to the magnetic model 

While magnetic models can simply explain bilateral SNRs which have a circular rim 
and two opposed limbs of similar shape and brightness (e.g. G003.8-00.3 and G327.6-I-14.6), 
we now argue that difficulties are encountered when one considers whether such models can 
also account for the detailed morphologies of other objects in Table ^. As will be discussed 
in Section [6.2.4[, an alternative explanation is the effect of the ISM. 
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6.2.1. Asymmetric SNRs 



There are objects in Table ^ which ahgn closely with the Plane, but for which the flux 
density and morphology of the bright opposed limbs is very different. The most striking 
examples of this class are G320.4-01.2 (Figure ^ and G308. 8-00.1, but this asymmetry is 
also present to a lesser extent in G078.2+02.1, G127.1+00.5, G350.0-02.0 (Figure |) and 
G356.3-01.5 (Figure |ll]). In the magnetic model, this morphology can be explained by 



a gradient in the field strength across the remnant ( Caswell fc Lerche 1979 ; [Reynolds fc 
Fulbright 19901) . 



Alternatively, for both G320.4-01.2 and G308. 8-00.1, it has been suggested that this 
appearance is due to the influence of a central pulsar on the surrounding shell ( [Manchester 



fc Durdin 1983[ ; [Kaspi et al. 1992t [Brazier fc Becker 1997|) , or because one side of the 
remnant encounters a region of higher density ([Seward et al. 1983| ; [Kaspi et al. 199"^ ). 
In the latter case, the low value of ip for these remnants suggests the presence of density 
gradients perpendicular to the Plane. There is no clear tendency for the limb closest to the 
Plane to be the brighter of the two. This is contrary to the suggestion of Caswell (1977), 
although he noted that this may only apply for remnants at large \z\. 



6.2.2. Distorted SNRs 

Some objects are signiflcantly distorted from a circular or elliptical shape, for example 
G350.0-02.0 (Figure g) and G166.0+04.3 (Figure 1). Observations of the latter at radio 
(Landecker et al. 1982, 1989), X-ray ( [Burrows fc Guo 1994| ; [Guo fc Burrows 1997[ ) and 



optical ( [Pineault et al. 1985| ) wavelengths, are all consistent with the model of Pineault 
et al. (1987). As discussed in Section [4.2[ , this explains the morphology in terms of 
the shock encountering a hot low density tunnel lying along the Galactic Plane, before 
propagating into dense material beyond. We have noted that G350. 0-02.0 has a similar 
appearance, which by analogy might also be explained in terms of such a model. The 
gross morphology of these remnants can be reproduced by magnetic models, but an abrupt 
change in density by a factor of 100, along an interface parallel to the Galactic Plane, is 
still required ( [Mineshige fc Shibata 1990[) . 



6.2.3. Elongated SNRs 



Some objects in Table ^ are signiflcantly extended along their symmetry axis, most 
notably G356.3-01.5 (Figure 0) and G296. 5+10.0 (Figure |^). This appearance can 
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be explained by tension associated with the ambient magnetic field, which would cause 
remnants to preferentially expand along field lines ( |lnsertis fc Rees 19911 |Mineshige! 



Shibata, fc Shapiro 1993| ; [Kozyczka fc 'I'enorio-'l'agle 1995| ). However for typical values of 



the explosion energy and ambient field strength, the shape of a SNR will not be appreciably 
distorted by the magnetic field until it reaches a diameter of ;^100 pc ([Mineshige fc Shibata| 



|1990| ; [Norman 1993| ), at which stage a remnant is nearing the end of its observable lifetime 



(e.g. [Shull, Fesen, fc Saken 19"89| ). In order for tension in the magnetic field to account for 
the observed extension, one must invoke a low explosion energy 10^° erg) and/or a high 
magnetic field 100 /iG). An alternative magnetic explanation for this elongation has 
been proposed by Roger et al. (1988). In the case of G296. 5+10.0, they speculate that the 
SNR shock may not couple to newly swept up material when the magnetic field is parallel 
to the direction of propagation. As the remnant evolves the shock speed is dependent on 
the angle with the magnetic field and the remnant will elongate. 

Bisnovatyi-Kogan et al. (1990) consider a non-magnetic explanation for G296. 5+10.0, 
wherein both the elongation and the bilateral appearance of the SNR are the result of it 
having expanded into a similarly extended pre-existing cavity. For this mechanism to apply 
for SNRs with low values of ■?/', however, such cavities must lie with their major axis along 
the Plane. 



6.2.4. An ISM model 

In the above discussion, we have considered three classes of objects which can be 
more comfortably accounted for by the surrounding ISM than by the magnetic model. 
We have considered objects with one side brighter than the other as might be due to a 
density gradient, a breakout morphology due to the presence of hot tunnels, and elongated 
remnants which might be due to an extended cavity. In order for any of these explanations 
to be valid, the observed distribution of ■0 requires that structures in the ISM must be 
consistently parallel to the Plane. 

As a mechanism for the bilateral appearance, we propose that the Galactic magnetic 
field produces density stratifications parallel to the Galactic Plane. Evidence to support 
this hypothesis is now presented. 

The first issue to consider is whether stellar wind-bubbles can be elongated by the 
surrounding field. Consider a homogeneous ISM of gas density nn, in which there is a 
uniform magnetic field of strength Bq. The bubble produced by an isotropic wind will 
expand more rapidly along field lines than across them, and will begin to be significantly 



- 18 - 



distorted when the internal pressure in the bubble falls below the external magnetic pressure 
( [Konigl 1982| ). Adapting the result of Stone & Norman (1992) for a wind from a proto-star, 
this will occur after a time ti, given by : 



X -5/2 / ^ \ 1/2 , , 3/4 



\Myr J \3 fiG J ylO^^ergy Vcm ^ 

where M is the mass of the star, Vw is the velocity of the wind and = is the 

associated mechanical luminosity. Therefore as a rough criterion, a wind bubble will become 
significantly elongated along ambient field lines if to ^ ^i, where to is the main-sequence 
lifetime of the star. 

Chevalier & Liang (1989) tabulate the properties of stellar wind-bubbles for different 
spectral types, from which we find that in a typical environment {Bq = 3 /iG, tlh = ^ cm~'^), 
a main-sequence wind-bubble will be dominated by the ambient field for a star of spectral 
type 09 or later. Such stars comprise the vast majority of the stellar population, and 
indeed Lozinskaya (1988) notes that wind-blown bubbles are often elongated. Since 
ti oc Bq^^'^ , only a slightly higher field than assumed above (^ 10 fiG) is required to distort 
wind-bubbles produced by even the most massive stars. Note that the extra energy input 
in the late phases of stellar evolution may complicate the simple picture presented here. 

On larger scales, Heiles (1979) has found that Hi shells within 10° of the Plane appear 
elongated in the longitudinal direction, which he suggests could also be due to the ambient 



field. Also, simulations of magnetised superbubbles (Tomisaka 1990, 1992; [Mineshige^ 
Shibata, fc Shapiro 1993|) and old SNRs ( [Norman 1993|) show that such structures also 
become elongated along field lines, but over larger length and time scales than for stellar 
wind bubbles. 

To summarise, we propose that where the magnetic field is oriented parallel to the 
Galactic disc, structures will be produced extended in the longitudinal direction, on 
approximately the same scale as SNRs. Thus close to the Plane, the structure of the 
Galactic ISM might take the form of tunnels, gradients and cavities all elongated in the 
disc of the Galaxy. When a SNR interacts with this structure, a bilateral appearance is 
produced (as in the model proposed by Bisnovatyi-Kogan, Lozinskaya, fc Silich 1990|) , with 



the axis aligned parallel to the Plane. 

A remnant formed from a massive progenitor will interact at its earliest stages with 
a complex and extensive GSM. During this period, a SNR can take on some aspects of 
bilaterality ([Igumenshchev, Tutokov, fc Shustov 199^ ; [Blondin, Lundqvist, fc Ghevaliei 



|1996| ; [Zhang, Wang, fc Chen 1996| ; [Gaensler et al. 1997[ ), but no alignment with the 
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Galactic Plane will be observed (unless the progenitor wind-bubble is elongated along the 
Plane through the process discussed above). Only once the shock propagates beyond its 
CSM and begins to interact with the ambient ISM should any trend in ifj manifest itself. 
Type la SNRs such as G327.6-I-14.6 may interact with the ambient medium at a much 
earlier stage, however. At the other end of the scale, as a remnant ages and expands, it 
becomes more and more likely to encounter complex inhomogeneities in the ISM which will 
distort its appearance. Hence, this mechanism for the bilateral appearance is most relevant 
for middle-aged remnants, rather than young or old ones. The uncertainty in SNR age 
estimates prevents us from demonstrating this property in the observed sample. However 
rough estimates for the ages of SNRs in Table ^ (using the E — t relation; Caswell & Lerche 
1979) give that most are between 1000 and 10 000 years old. 



6.3. Exceptional cases: SNRs G296. 5+10.0 and G327.6+14.6 

SNRs G296.5+10.0 (Figure 0) and G327.6+14.6 (Figure |) are discrepant in their 
values of both being oriented at almost right angles to the Plane. However, these two 
remnants are also notable in that they are both at a significant height above the Plane 
(175-350 pc and 450 pc, respectively; [Roger et al. 1988) ). 

For bilateral SNRs at low latitudes, we have proposed ISM structures running parallel 
to the Plane. At higher latitudes, the ISM is dominated by perpendicular structures in 
the form of chimneys ( [Norman fc Ikeuchi 1989| ; [Normandeau, Taylor, fc Dewdney 1996[ ) 



and "worms" ([Heiles 1984[ ; [Koo, Heiles, fc Reach 1992]) . Stellar polarization measurements 



demonstrate the existence of large-scale loops of magnetic field also extending perpendicular 
to the Plane ( [Mathewson fc Ford 1970[ ). Thus that bilateral SNRs might tend to have 
ip ~ 90° at large distances from the Plane is consistent with the extrinsic models described 
in Sections ^ and ^:T^ SNRs G296.5+10.0 and G327.6+14.6 can therefore still be 
incorporated within the proposed model. One can speculate that bilateral SNRs at 
intermediate scale heights might be in transition between remnants with ~ 0° and those 
for which ifj ~ 90°. However the uncertainty in the distances to SNRs in our sample 
prevents us from making meaningful conclusions on any latitude dependence at present. 

These two SNRs are regarded as the prototypes for the bilateral morphology (e.g. 
KC87; [Roger et al. 19S5[). That the most striking examples of bilaterality constitute 



50% of known remnants for which |6| > 10° may be indicative of the homogeneity of the 
environment far from the Plane. It is highly likely that the radio all-sky surveys currently 
being carried out ([Large et al. 199*4 ; CJondon et al. 19"96 ) will discover new high latitude 



SNRs. Whether such remnants also have a well-ordered appearance and a high degree 
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of symmetry will be useful (but not conclusive) evidence regarding the model we have 
proposed. 

We note that intrinsic models have also been proposed for these two SNRs: Storey 
et al. (1992) argue that the symmetry present in G296. 5+10.0 over a wide range of spatial 
scales generally argues against magnetic models for this object, while Willingale et al. 
(1996) hypothesise that a central compact object with opposed jets could generate the 
morphology of G327.6+14.6. Whether these remnants represent an extrinsic mechanism 
operating at high Galactic latitude, or intrinsic effects, either way, their high value of if) is 
not at odds with the alignment observed in the remainder of the sample. 



6.4. Explanations for non-bilateral structure 



As noted in Section f).2\ the vast majority of remnants do not have a clear bilateral 
appearance. The wide range of SNR morphologies observed (clearly demonstrated in 
the catalogue of WG96) is indicative that complex, small-scale effects can dominate the 
appearance of SNRs. Specifically, the absence of bilaterality can be explained by any of the 
following conditions: 



A dominant component of the field along the line of sight, so that the projection of 
the emitting region is no longer two distinct arcs of emission. KC87 estimate that a 
SNR with a cylindrical 3D geometry will appear as a uniform ring of emission with no 
bilateral appearance ii 6 < 25°, and will be "confused" if 25° < 9 < 45°. This results 
in 24% (rounded to 30% in KC87) of SNRs having a non-bilateral appearance, simply 
due to orientation. 

An inhomogeneous ISM. In this case, encounters with molecular clouds, cavities, and 
other inhomogeneities distort a remnant sufficiently so that no symmetrical structure 
can be observed (e.g. [Pineault et al. 1993| ; pieynolds fc Moffett 1993| ; [Koo fc Moon 



19971) . 



A disordered ambient magnetic field, so that in the magnetic model, no position angle 
around the circumference has enhanced emissivity. It is suggestive that some bilaterals 
are highly polarized (e.g. [Lalitha et al. 195^ ; [Milne fc Haynes 1994) ), indicating a 
well-ordered magnetic field. 
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7. Conclusion 

We have presented 1.4 GHz radio observations of two southern supernova remnants, 
G003.8-00.3 and G350.0-02.0. The morphologies of these two SNRs are distinct — 
G003.8-00.3 shows a filamentary two-arc appearance in a classic "barrel" shape, while 
G350.0-02.0 strongly resembles the breakout morphology of G166.0+04.3 (VRO 42.05.01). 

Despite their differences, both remnants have a clear bilateral appearance. Furthermore, 
the axis defined by this classification in both cases aligns almost exactly with the Galactic 
Plane. On examination of the entire sample of known Galactic SNRs, we find 17 SNRs with 
a high level of bilateral morphology, of which nine are oriented to the Plane at 12° or less. 
With a minimum of assumptions regarding the three-dimensional geometry, the probability 
of this occurring by chance is 7 x 10~^. Considering SNRs G296. 5+10.0 and G327.6+14.6 
as distinct from the remainder of the sample reduces this probability even further. 

This trend rules out "intrinsic" models for the bilateral morphology, all of which involve 
an axis defined by the orientation or motion of the progenitor star. Such models are not 
expected to produce remnants which align in any consistent way with the Galactic Plane. 

The general alignment of bilateral remnants with the Plane can be explained by 
"magnetic" models, whereby an ambient magnetic field generates bright regions around the 
SNR shell where the shock normal is perpendicular to the field. This is a result of either 
enhanced field compression or preferential electron acceleration in these regions. However, 
distorted morphologies in bilateral SNRs can be better explained if the ISM into which 
these remnants expand consists of structures preferentially elongated in the longitudinal 
direction. We propose that this occurs in the form of expanding wind-bubbles distorted by 
the surrounding magnetic field. 

The majority of SNRs do not have a distinct bilateral appearance. For some SNRs, 
this can be explained as an orientation effect, when the bilateral axis is close to the line of 
sight. The effects of a disordered magnetic field and an inhomogeneous ISM can further 
account for the observed diversity in morphologies. 

This study has demonstrated the value of high resolution observations of Galactic 
SNRs. Further observations, in particular of the polarization and environment of bilateral 
SNRs, are needed to test the proposed models. Undoubtedly there also remain more 
examples of bilateral SNRs to be discovered, particularly in the first quadrant and at high 
Galactic latitudes, which will improve the available statistics. 



While theoretical modelling of the evolution of stellar wind bubbles is quite sophisticated 
( [Garcia- Segura, Langer, fc Mac Low 1996| and references therein), the effect of external 
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magnetic fields is generally ignored. The simple calculations made here suggest that such 
fields can have a significant effect on the dynamics of a bubble and should be included in 
subsequent modelling. Similarly, the interaction of SNRs with surrounding material usually 
assumes spherical symmetry ( |(Jhevalier fc Liang 1989| ; I'enorio-'l'agle et al. 1991 ); detailed 



simulations of the interaction of remnants with the (non-spherical) structures proposed here 
would be highly desirable. 

An interesting extension of the results described would be to determine whether 
bilateral SNRs in other edge-on spiral galaxies are aligned along preferred position angles, 
and whether these angles correlate with measurements of the corresponding galactic 
magnetic field. Of particular interest is NGC 4631 which largely has a magnetic field 
oriented perpendicular to the disc ( [Hummel, Beck, fc Dahlem 1991| ; [Golla fc Hummel 1994| ). 
At present, existing observations of SNRs in other spiral galaxies (e.g. puric et al. 1995| ; 



IMuxlow et al. 1994|) are of insufficient resolution to derive a bilateral subset. However, the 



planned upgrade to the VLA ([Bastian fc Bridle 1996|) will provide sufficient improvement 
in sensitivity and resolution to produce images of SNRs in these galaxies of a quality 
comparable to existing Galactic observations. 
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FIGURE CAPTIONS 

Fig. 1.— 1.4 GHz VLA image of SNR G003.8-00.3. The linear greyscale is from -0.5 to 
3.0 mJy beam~^, and the peak surface brightness in the image is 15.8 mJy beam~^. The 
synthesised beam, shown at the lower right of the Figure, corresponds to a resolution of 
15" X 9". The solid line represents a constant Galactic latitude b = —0°A. The dashed line 
represents the axis of symmetry for this SNR, as discussed in the text. 

Fig. 2. — Image of SNR G350. 0-02.0, produced by combining 1.4 GHz data from the VLA 
and from Parkes. The linear greyscale is from -2.4 to 7.5 mJy beam~^, and the peak surface 
brightness in the image is 58.5 mJy beam~^. The synthesised beam, shown at the lower right 
of the Figure, corresponds to a resolution of 21" x 18". A line of constant Galactic latitude 
and the axis of symmetry for the SNR are shown as in Figure |1|. Regions A, B and C, as 
discussed in the text, are indicated. 

Fig. 3.— 1.4 GHz image of SNR G166.0+04.3, reproduced from Pineauh et al. (1987) 
and precessed to J2000 coordinates. A line of constant Galactic latitude and the axis of 
symmetry for the SNR are shown as in Figure ||. The "wing" and "shell" components as 
defined by Landecker et al. (1982) are shown. 
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Fig. 4.— 0.843 GHz image of SNR G315.4-02.3 (ROW 86), reproduced from WG96. 

Fig. 5.— 0.843 GHz image of SNR G332.4-00.4 (RCW 103), reproduced from WG96. 

Fig. 6. — 1.4 GHz image of G046.8-00.3, reproduced from Dubner et al. (1996) and 
precessed to J2000. Annotations to the image are as in Figure |I|. 

Fig. 7. — 0.843 GHz image of SNR G296. 5+10.0, reproduced from KC87 and precessed to 
J2000. Annotations to the image are as in Figure |1|. 

Fig. 8.— 0.843 GHz image of SNR G320.4-01.2, reproduced from WG96. Annotations to 
the image are as in Figure [^. 

Fig. 9. — 0.843 GHz image of SNR G327.6+14.6, reproduced from KC87 and precessed to 
J2000. Annotations to the image are as in Figure |l|. 

Fig. 10.— 0.843 GHz image of SNR G332.0+00.2, reproduced from WG96. Annotations to 
the image are as in Figure 

Fig. 11.— 0.843 GHz image of SNR G356.3-01.5, reproduced from Gray (1994). 
Annotations to the image are as in Figure |^. 

Fig. 12. — Histogram of the distribution of ■0, the acute angle between the bilateral axis of 
a SNR and the Galactic Plane, for SNRs listed in Table 

Fig. 13. — Representation of a bilateral SNR in Cartesian space. The vector corresponds 
to the axis of the remnant, assuming cylindrical symmetry. The observed angle between the 
projected axis and the Galactic Plane is -0 = \(j)\. 
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Date 


Array 


Source 


Time on source 


1996 Jan 20 


CnB 


G003. 7-00.2 


165"^ 






G350.0-01.8 


90"^ 


1996 May 21 


DnC 


G350.0-01.8 


24m 


1996 Jun 18 


DnC 


G003.7-00.2 


31™ 






G350.0-01.8 


24m 


1996 Jun 22 


DnC 


G003.7-00.2 


47m 






G350.0-01.8 


40™ 



Table 1: VLA observations. 
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Observing date 


1996 Sep 07 


Centre frequency 


1.40 GHz 


Bandwidth 


40 MHz 


System temperature 


25 K 


RMS noise 


24 mJy beam"-^ 


Beam width (FWHM) 


14^85 ± 0^25 


Map gridding 


6' 



Table 2: Parkes observations. 
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G003.8-00.3 


G350.0-02.0 


VLA pointing centres^ 


17''55™30^ -25°51'00" 


(1) 17''28™30^ -38°37'30" 






(2) 17''27™20^ -38°27'30" 






(3) 17''26™55^ -38°10'00" 






(4) 17^26^00% -38°23'00" 


Resolution 


15" X 9", PA = 60° 


21" X 18", PA = -13° 


RMS noise (mJy beam~^) 


0.06 


0.15 


Derived centre of SNR^ 


17'^55'"28^ -25°50' 


17^'27'"50^ -38°30' 


Derived centre of SNR^ 


003.78, -00.29 


349.95, -02.03 


Diameter 


13' 


47' X 44' 


1.4 GHz flux density (Jy) 


1.7 ±0.1 


22.3 ±0.3 


Spectral index 


-0.65 ± 0.05 


-0.4 ±0.1 


'S'lGHz (Jy) 


2.3±0.1 


25.7±0.9 


Si GHz (10-2^ W m-2 sr-i) 


2.0±0.1 


1.7±0.1 



Table 3: Observational, measured and derived parameters for each SNR. 



^Equatorial coordinates a, 5 (J2000) 
^Galactic coordinates I, b 
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SNR 


Other Name 




Figure 


Reference 


G003.8-00.3 




6° ± 2° 


1 


This paper 


G046.8-00.3 




21° ±3° 


1 


1 


G078.2+02.1 


7 Cygni 


10° ± 2° 




2 


G093.3+06.9 


DA 530 


15° ± r 




3 


G127.1+00.5 




3° ±3° 




4 


G156.2+05.7 




64° ± 4° 




5 


G166.0+04.3 


VRO 42.05.01 


2° ±2° 


1 


6,7 


G296.5+10.0 


PKS 1209-51/52 


80° ± 3° 





8, 9 


G302.3+00.7 




42° ± 2° 




10 


G308.8-00.1 




9° ±3° 




10 


G318.2+00.1 




12° ±3° 




10 


G320.4-01.2 


ROW 89, MSH 15-5^ 


7° ±3° 


1 


10 


G327.6+14.6 


SN 1006 


86° ± 2° 


1 


8, 9 


G332.0+00.2 




25° ± 4° 





10 


G336.7+00.5 




46° ± r 




10 


G350.0-02.0 




5° ± 2° 





This paper 


G356.3-01.5 




2° ±2° 


11 


11 



Table 4: Bilateral SNRs and their orientation with respect to the Galactic Plane. 



References. — (1) Dubner et al. (1996) (2) Pineault & Chastenay (1990) (3) Lalitha et al. (1984) 
(4) Joncas, Roger & Dewdney (1989) (5) Reich, Fiirst and Arnal (1992) (6) Landecker et al. (1982) 
(7) Pineault, Landecker & Routledge (1987) (8) KC87 (9) Roger et al. (1988) (10) WG96 (11) Gray (1994) 
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